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We present magic wavelengths for the nS - nPi/ 2 ^ 3/2 and nS - mD 3 / 2 ,s /2 transitions, with the 
respective ground and first excited D states principal quantum numbers n and m, in the Mg'^, Ca’'", 
Sr+ and Ba'*’ alkaline earth ions for linearly polarized lights by plotting dynamic polarizatbilities of 
the nS, nPi/2,3/2 and 01,1)3/2,5/2 states of the ions. These dynamic polarizabilities are evaluated by 
employing a relativistic all-order perturbative method and their accuracies are ratified by comparing 
their static values with the available high precision experimental or other theoretical results. More¬ 
over, some of the magic wavelengths identified by us in Ca^ concurs with the recent measurements 
reported in [Phys. Rev. Lett. 114,223001 (2015)]. Knowledge of these magic wavelengths are 
propitious to carry out many proposed high precision measurements trapping the above ions in the 
electric fields with the corresponding frequencies. 

PACS numbers: 32.10 Dk,31.15.Dv, 31.15 ap 


State-insensitive trapping techniques have lead to 
tremendous advancements in the manipulation and con¬ 
trol of atoms in far detuned optical traps. In this ap¬ 
proach, the atoms are trapped at the wavelengths (re¬ 
lated to frequencies) of an external electric field at which 
the differential light shift of an atomic transition, that 
is intended to be probed, due to the Stark effects nul¬ 
lify. These wavelengths are specially referred to as magic 
wavelengths (AmagicS M- It has been demonstrated ear¬ 
lier ability of trapping neutral atoms inside high-Q cav¬ 
ities at AmagicS in the strong coupling regime, which is 
important in the quantum computation and communica¬ 
tion schemes [3]. This technique is now widely used to 
carry out many high precision measurements by eliminat¬ 
ing large systematics due to stray electric fields. Another 
notable application of these wavelengths is to perform 
clock frequency measurements Q, especially for optical 
frequency standards [^, that are in turn useful to probe 
temporal and spatial variations of the fundamental con¬ 
stants @ and improving global positioning systems 
Knowing AmagicS of atomic systems are also very useful in 
the field of quantum state engineering [7(|, extracting out 
precise values of the oscillator strengths Q, etc. because 
of which extensive studies, both experimentally and the- 
oretically, have been carried out in many atoms recently 
d, l9l-[T^. On the other hand, singly charged alkaline 
earth ions are advantageous to carry out very high preci¬ 
sion measurements using ion trapping and laser cooling 
techniques. Some of the prominent examples are the op¬ 
tical frequency standards HP, probing variation of 
fundamental constants [3, Q i parity non-conservation ef¬ 
fects [H, [1^ etc. Advantages of these ions owe to their 
metastable D states that provide longer probe times dur¬ 
ing interrogation of measurements. To reduce system- 
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atics in these measurements, state insensitive measure¬ 
ments can be more pertinent that will require knowl¬ 
edge of AmagicS in these ions. In fact, AmagicS for the 
45' —5> 31)3/2,5/2 and 45 —>■ 4Pi/2,3/2 transitions in Ca'*' 
are recently observed m. 

In this Rapid Communication, we report AmagicS for 
the nS — nPii 2 ^ 3 i 2 and nS — transitions in the 

Mg+j Ca“'", Sr"*" and Ba"*" ions, for the principal quantum 
numbers of the ground state n and of the first excited D 
states m of the respective ions, for the commonly used 
linearly polarized lights in the experiments by evaluating 
dynamic dipole polarizabilities of these ions accurately. 

Dominant contribution to the change in energy of a 
state with angular momentum and its com¬ 

ponent m„ of an atomic system due to interaction with 
an external electric field 5(w) of frequency u is given by 

^Pstark ~ (^) 

where a„(a;) is the dynamic dipole polarizability of the 
atomic state jjnjWn) and given by 


a„(w) = ai^^w) -b 


3to| - j„(j„-b 1) ^2) 

Jn(2j„ - 1) 


(w), 


( 2 ) 


for the scalar and tensor components as al°^(a;) and 
(w), respectively. Expressions for these quantities in 
terms of the reduced electric dipole (El) matrix elements 
can be found in Refs. iSIl- As discussed in our previous 
work [1^, we estimate al^^s (for A: = 0,2) by expressing 

+ al^'> {v) (3) 


with the notations c, vc and v representing the contribu¬ 
tions due to the core, core-valence and valence correlation 
effects, respectively. Dominant valence correlations are 
estimated by calculating electric dipole (El) matrix el¬ 
ements between many low-lying intermediate states and 
referred to as “Main” result. This is evaluated using 
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FIG. 1: (Color online) Dynamic polarizabilities (in an) of the 
351/2, 3 Pi/2 , 3/2 and 31 ) 3 / 2 , 5/2 states in Mg’^. Arrows indicate 
positions of magic wavelengths (Amagics) in nm. 


an all order relativistic coupled-cluster method with the 
singles and doubles approximation (SD method) as de¬ 
scribed in [iilig. Other small contributions from the 
higher excited states (referred to as “Tail” contribution), 
core and core-valence correlations are estimated using 
lower order perturbative methods as described in [l8| . To 
construct the single particle orbitals for the SD method, 
we have used total 70 B-spline functions with a cavity of 
radius R = 220 au. However, a few El matrix elements 
involving the states in Ba+ are taken from [2l| 

to use as many as El matrix elements for more accurate 
evaluations of the Main contributions to a„s of the 5Z) 
states of this ion. Estimated static polarizabilities in the 
above procedure of the considered ground and excited 
states of the alkaline earth metal ions using these matrix 
elements and experimental energies from national insti¬ 
tute of standards and technology (NIST) [ 2 ^ are listed 
in the Supplemental material along with the individual 
core, core-valence and Tail contributions. These values 
are also compared with the available precise experimen¬ 
tal results and theoretical calculations. We find excellent 
agreement between our results with these values imply¬ 
ing that our calculations are precise enough to predict 
the magic wavelengths in the considered ions reliably. 

In Tables U [III [TIT] and |IVl we list AmagicS for the 
nS — nPii 2 ^ 3 i 2 and nS — transitions in the 

Mg+, Ca’*', Sr+ and Ba’*' ions respectively. They are ob¬ 
tained by locating the crossings between the dynamic po¬ 
larizabilities of the nS, nPi/ 2 , rLP^/ 2 , rnD^/ 2 , and mD ^/2 
states of Mg+,Ca''", Sr+ and Ba"*" plotted against the fre¬ 
quencies as shown in Figs. [Tl[51[3]and|3]respectively. The 
arrows in the figures indicate the positions of Amagic(avg) 
values which are the average of the magic wavelengths 
at different rrij sub-levels. The Aniagic(avg) are not pre¬ 
sented for cases where magic wavelengths are (a) not 
identified for all rrij sub-levels, (b) supporting different 
kind of trapps for all mj sublevels (c) separated by large 


TABLE I: Magic wavelengths (Amagics) in nm and correspond¬ 
ing polarizabilities (in au) in the Mg^ ion. List of the resonant 
wavelengths (Aress) and AmagicS are also given sequentially in 
increasing order to indicate their respective placements. 


Resonance 

-^res 

\mj\ 

'^magic 

<^(Ainagic) 

Transition 


35 

“ 3Pi/2 


37i/2“37)3/2 

279.0777 

1/2 

276.12 

-1267.54 

3Fi/2-35i/2 

280.2704 




3Fi/2-45i/2 

292.8683 

1/2 

307.23 

202.59 

Transition 


35 

“ 3 P 3/2 


3 B 3/2 “ 3D3/2 

279.7930 




37*3/2 ~ 3D5/2 

279.7998 




37*3/2 ~ 35i/2 

279.5528 




37*3/2 — 45i/2 

293.6510 

1/2 

323.94 

136.51 

Transition 


35 

— 37)3/2 


37 ) 3/2 — 57^5/2 

310.4809 

1/2 

313.73 

169.42 



3/2 

314.18 

167.53 

37 ) 3/2 — 5 P 3/2 

384.8335 

3/2 

385.01 

72.56 

37 ) 3/2 — 5Pi/2 

385.0385 

1/2 

385.38 

73.47 

37 ) 3/2 — 47^5/2 

448.1327 

1/2 

705.03 

41.44 



3/2 

900.91 

38.67 

37)3/2 — 4P3/2 

1091.5270 




37 ) 3/2 — 4Pi/2 

1095.1779 




Transition 


35 

~ 37)5/2 


37 ) 5/2 — 57^7/2 

310.4715 




37 ) 5/2 — 57^5/2 

310.4721 

1/2 

313.67 

169.66 



3/2 

313.80 

169.01 



5/2 

314.39 

166.70 

37 ) 5/2 — 5 P 3/2 

384.8209 

3/2 

385.12 

73.58 



1/2 

385.17 

73.56 

37 ) 5/2 — 47^7/2 

448.1130 




37 ) 5/2 — 47^5/2 

448.1150 

1/2 

699.0 

41.57 



3/2 

737.27 

40.79 

37)5/2 — 4P3/2 

1091.423 






wavelength(nm) 
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FIG. 2: (Color online) Dynamic polarizabilities (in au) of the 
45 i/2, 4Pi/ 2,3/2 and 31 ) 3 / 2 , 5/2 states in Ca"*". Locations of 
average AmagicS (in nm) are indicated by arrows. 


wavelengths for different rrij sub-levels. We also give the 
total an values at the corresponding AmagicS considering 
different mj values in the tables. Before we proceed in 
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TABLE II: AmagicS (in nm) and corresponding polarizabilities 
(in au) in the Ca"^ ion. Ares are also presented to identify 
the placements of AmagicS clearly between two resonant wave¬ 
lengths. Results from this work ([T]) are compared with the 
findings from Ref. @]. 


Resonances 

Ares 

\mj\ 

Amagic 

Q:( Amagic) 

Reference 

Transition 



45- 

4Pi/2 


4Pi/2 

— 4I?3/2 

315.887 

1/2 

368.10 

-485.18 

[T] 





368.0149 

-477.2554 

[8] 

4Pi/2 

- 55i/2 

370.603 

1/2 

395.18 

2933.47 

[T] 



395.1807 

2896.7954 

[8] 

4Pi/2 

— 45i/2 

396.847 

1/2 

697.65 

110.39 

[T] 





690.817 

110.0849 


to 

~ 3I?3/2 

866.214 




Transition 



45- 

4 P 3/2 


4 P 3/2 

~ 4I?5/2 

317.933 





4 P 3/2 

— 4I?3/2 

318.128 

1/2 

369.72 

-520.85 

[T] 





369.6393 

-512.1655 

[8] 

4 P 3/2 

- 55i/2 

373.690 




4 P 3/2 

— 45i/2 

393.366 

3/2 

395.77 

125.27 

[T] 





395.774 

123.4906 

[8] 




1/2 

396.23 

-3118.47 

[T] 





396.2315 

-3077.3881 

[8] 

45i/2 

— 4Pi/2 

396.847 

3/2 

678.35 

113.38 

[T] 





672.508 

113.0150 





1/2 

693.76 

110.96 

[T] 





687.022 

110.6606 


4 P 3/2 

~ 3I?3/2 

849.802 

1/2 

850.12 

81.57 

[T] 




3/2 

850.92 

96.09 

[T] 

4 P 3/2 

~ 3I?5/2 

854.209 





Transition 



45- 

3 P 3/2 


45i/2 

“ 4 P 3/2 

393.366 

1/2 

395.80 

3.13 

[T] 





395.7981 

1.1711 





3/2 

395.79 

8.02 

[T] 





395.7970 

7.1019 


45i/2 

— 4Pi/2 

396.847 





37?3/2 

~ 4 P 3/2 

849.802 

1/2 

850.33 

96.12 

[T] 





850.335 

95.0011 

[8j 

37?3/2 

— 4Pi/2 

866.214 

3/2 

884.54 

94.23 

[T] 





887.382 

92.9908 

[8] 




1/2 

1252.44 

84.17 

[T] 





1308.590 

82.4644 

[8j 

Transition 



45- 

3 P 5/2 


45i/2 

— 4 P 3/2 

393.366 

1/2 

395.79 

2.44 

[T] 





395.7982 

0.5371 





3/2 

395.79 

3.49 

[T] 





395.7978 

3.1792 





5/2 

395.79 

8.83 

[T] 





395.7968 

8.4633 


45i/2 

— 4Pi/2 

396.847 




37?5/2 

~ 4 P 3/2 

854.209 

3/2 

1052.26 

88.06 

[T] 





1074.336 

86.4837 





1/2 

1271.92 

83.90 

[T] 





1338.474 

82.1167 

[8] 


presenting AmagicS, we would like to clarify that a^s of the 
nS states can have two important resonant transitions 
{nS — nPij 2 and nS — nPs/ 2 ) in the wavelength range 
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FIG. 3: (Color online) Dynamic polarizabilities (in au) of the 
55 'i/2, 5Pi/ 2,3/2 and 41 ) 3 / 2 , 5/2 states in Sr"^. Arrows indicate 
positions of the average Amagic in nm. 


considered in this work. Therefore, anS of the nS states 
are generally small except in the close vicinity of these 
transitions. In contrast, a-nS of the nP and mD states 
can have significant contributions from several resonant 
transitions in the considered wavelength range. Thus, 
they are expected to cross with of the nS states in 
between these resonant transitions. We also list the res¬ 
onant wavelengths (Ares) in the above tables to highlight 
respective placements of AmagicS- We find similar trends 
in locating magic wavelengths for all the considered ions, 
except for few cases where AmagicS are missing. 

Mg+: We are able to locate two AmagicS magic wave¬ 
lengths for the 35 — 3Pi/2 transition in Mg+, one just 
before the 3Pi/2 — 31? resonant transition at 276.12 nm 
and one just after the 3Pi/2 — 45 resonant transition 
at 307.23 nm. The first Amagic supports a blue detuned 
trap, whereas the second one a red detuned trap. How¬ 
ever, only one Amagic for the 35 — ^P 3 / 2 {'m-j = il/2) 
transition is located at 323.94 nm and nothing is found 
for the 35 — 3P3/2{mj = ±3/2) transition. On the other 
hand, three AmagicS are identified for the 3I?3/2 —35 tran¬ 
sition placed between six different resonances. First one 
is located around 314 nm between the 3I?3/2 — 555/2 and 
37?3/2 ~ 5P3/2 resonant transitions. The next Amagic is 
located at the sharp intersection of polarizability curves 
close to the 3Z?3/2 — 5P3/2 and the 3Z?3/2 — 5Pi/2 reso¬ 
nances as seen in Fig. [TJ The last Amagic is located at 
705.03 nm for rrij = ±1/2 and 900.91 nm for rrij = ±3/2 
sub-levels, hence it is of limited experimental use. For 
the 35 — 3I?5 /2 transition, only one Amagic appears for all 
rrij sub-levels around 314 nm. All the observed AmagicS 
along with the AresS are tabulated in Table [H It can be 
noticed from this table that most of these magic wave¬ 
lengths favor weak red detuned traps. 

Ca+: At least three AmagicS are located for the 4Pi/2 — 
45 transition in Ca"*" at 368.10, 395.18 and 697.65 nms 
between the 4Pi/2 — 4Z?3/2, 4Pi/2 — 55, 4Pi/2 — 45 and 
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TABLE III: AmagicS (in nm) and their corresponding polariz¬ 
abilities (in an) in Sr"*' and Ares contributing to the polariz¬ 
abilities of the transitional states are also listed. 


Resonances 

4^res 

I'm A 

-^magic 

<^(Axnagic) 

Transition 


55 

- 5Pi/2 


5Pi/2 — 5 D 3/2 

338.0711 




5Pi/2 — 65i/2 

416.1796 




5Pi/2 — 55i/2 

421.5524 

1/2 

769.44 

126.34 

5Pi/2 — 4 D 3/2 

1091.4874 




Transition 


55 

— 5 P 3/2 


5 P 3/2 — 5 D 5/2 

346.4457 

3/2 

347.38 

-204.86 

5 P 3/2 — 5 D 3/2 

347.4887 

1/2 

347.57 

-205.66 

5 P 3/2 — 5Si/2 

407.7714 

3/2 

416.62 

168.54 



1/2 

419.30 

-1521.68 

55i/2 - 5Pi/2 

421.5524 




5 P 3/2 — 65i/2 

430.5447 

1/2 

438.37 

817.75 



1/2 

716.72 

134.25 



3/2 

724.92 

132.83 

5 P 3/2 — 4 D 3/2 

1003.6654 

1/2 

1004.47 

108.08 



3/2 

1009.80 

108.90 

5 P 3/2 — 4 D 5/2 

1032.7309 




Transition 


55- 

- 47 ) 3/2 


55i/2 - 5 P 3/2 

407.7714 

3/2 

417.00 

17.36 



1/2 

417.00 

13.83 

55i/2 — 5Pi/2 

421.5524 




47 ) 3/2 — 5 P 3/2 

1003.6654 

1/2 

1005.30 

108.93 



3/2 

1082.38 

106.29 

47 ) 3/2 — 5Pi/2 

1091.4874 




Transition 


55- 

“ 47 ) 5/2 


55i/2 - 5 P 3/2 

407.7714 

1/2 

417.01 

12.16 



5/2 

417.00 

18.43 



3/2 

417.00 

14.25 

55i/2 — 5Pi/2 

421.5524 




47 ) 5/2 — 5 P 3/2 

1032.7309 
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FIG. 4: (Color online) Dynamic polarizabilities (in an) of 
the 6 S'i/ 2 , 6 Pi /2 53/2 ^-nd 51 ) 3 / 2 , 5/2 states in Ba"*" with their 
average AmagicS (in nm) pointed by arrows. 


4 ^’i /2 ~ 31 ) 3/2 resonant transitions. A good agreement 
is found between our results and findings by Tang et. 
al. @ except for our Amagic at 697.65 nm, which is red 


TABLE IV: AmagicS (in nm) and their corresponding polariz¬ 
abilities (in an) in Ba"*". AresS are also given to identify the 
placements of AmagicS between them. 


Resonances 

-^res 

\mA 

'^magic 

Qf('^magic) 

Transition 


65 

— 6 P 1/2 


6 P 1/2 — 85i/2 

264.72608 

1/2 

323.33 

-93.74 

6 P 1/2 — 67 ) 3/2 

389.17790 

1/2 

451.73 

-4574.49 

6 P 1/2 — 75 i /2 

452.4926 

1/2 

468.90 

953.80 

6 P 1/2 — 65i/2 

493.4077 

1/2 

599.28 

311.58 

6 P 1/2 — 57 ) 3/2 

649.6898 




Transition 


65 

— 6 P 3/2 


6 P 3/2 — 85i/2 

277.13528 

3/2 

348.85 

-133.47 



1/2 

360.63 

-158.70 

6 P 3/2 — 6 P 5/2 

413.06491 




6 P 3/2 — 6 P 3/2 

416.60014 

1/2 

416.66 

-460.38 



3/2 

416.06 

-459.32 

6 P 3/2 — 65i/2 

455.4033 

3/2 

475.72 

374.67 

6 P 3/2 — 75 i /2 

489.9927 

1/2 

552.60 

441.01 



3/2 

561.01 

406.42 

65i/2 — 6 P 1/2 

493.4077 




6 P 3/2 — 5 P 3/2 

585.3675 

1/2 

586.24 

336.10 



3/2 

593.40 

321.95 

6 P 3/2 — 5 P 5/2 

614.1713 




Transition 


65 

— 57 ) 3/2 


65i/2 — 6 P 3/2 

455.4033 

3/2 

480.44 

15.88 



1/2 

480.81 

-15.82 

65i/2 — 6 P 1/2 

493.4077 




57 ) 3/2 — 6 P 3/2 

585.3675 

1/2 

585.98 

336.18 



3/2 

592.46 

323.72 

57 ) 3/2 — 6 P 1/2 

649.6898 

1/2 

767.81 

194.05 

Transition 


65 

— 57 ) 5/2 


65i/2 - 6 P 3/2 

455.4033 

5/2 

480.26 

31.04 



3/2 

480.71 

-7.19 



1/2 

480.93 

-26.46 

65i/2 — 6 P 1/2 

493.4077 




57 ) 5/2 — 6 P 3/2 

614.1713 

3/2 

666.64 

239.13 



1/2 

718.18 

211.13 


detuned as compared to their value. Similarly, a total 
of four AmagicS are located for the 45” — 4 P 3/2 transition 
as shown in Fig. [5] and listed in Table m We also find 
AmagicS at 395.77 nm and 396.23 nm for the 45 — 4 P 3/2 
transition within the fine structure splittings of the 4P 
states which coincides with the experimental observa¬ 
tions at 395.7992(7) nm and 395.7990(7) nm respectively 
Similar to the 4Pi/2 — 45 transition, the other ob¬ 
served AmagicS at 678.35 nm and 693.76 nm are red de¬ 
tuned as compared to findings in Ref. 0. Moreover, they 
have missed a Amagic around 850 nm which lies in between 
the 4 ^ 3/2 — 31 ) 3/2 and 4 P 3/2 — 31 ) 5/2 resonant transitions. 
AmagicS for the 45 — 31 ) 3/2 and 45 — 31)5/2 transitions, de¬ 
tected between the fine structure splitting of the 4P state 
at 395.79 nm, offer a scope to trap Ca"*" at very small po¬ 
larizabilities and are not much of practical use. As seen in 
Fig. [2 Q!„s of the 31 ) 5 / 2 (m^- = ±5/2) state attain a con¬ 
stant value through out the wavelength range considered 
in this work. Consequently, only one crossing is noticed 
for this state at 395 nm when of the 45 state changes 
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sign close to the 45” — 4P resonances. The other identi¬ 
fied AmagicS at 1052.26, 1252.44 and 1271.92 nms for the 
4S' — 3 ZI 3/2 and 45 — 3115/2 transitions with different rrij 
values are blue detuned compared to the earlier findings 
of Ref. Q . The possible reasons for these disagreements 
could be that the plots for a„s of the transitional states 
cross at very small angles at these wavelengths thus mak¬ 
ing them very in-distinctive. Moreover, Tang et. al. have 
performed semi-empirical calculations of polarizabilities 
in the relativistic framework that do not take into ac¬ 
count electron correlation effects rigorously like our SD 
method. 

Sr+: From Fig. [Sj we are able to locate only one 
magic wavelength at 769.44 nm for the 55 - 5Pi/2 tran¬ 
sition of Sr+. As given in Table lllli AmagicS for the 
55—5P3/2 transition are, however, systematically placed 
between various resonances. We also locate a crossing 
between the 55 and = ±1/2) states at 438.37 

nm. All identified AmagicS for the 55 — 5P3/2 transi¬ 
tion are typically small in magnitude with exception for 
Amagic = 419.302 nm. At this Amagic, it is recommended 
to use blue tuned trap for the to/ = ±1/2 sub-levels of 
Sr+. For the 4P>3/2 —>■ 55 transition, Amagic = 1007.22 
nm lies in the infrared region and is recommended for red 
detuned trap. In Fig. [3l we also observe that a„s for the 
transitional states at Amagic = 417 nm for the 55 — 4113/2 
and 55 — 4ZI5/2 transitions are very small (almost ap¬ 
proach to zero). Therefore, it may not be advantageous 
to trap the Sr+ ion at this wavelength. 

Ba+: Four AmagicS are found for the 65 — 6 P 1/2 tran¬ 


sition in Ba+ at 323.33, 451.73, 468.90 and 599.28 nms 
among which Amagic = 451.73 nm is in direct vicinity 
of the 6P1/2 — 75 i /2 resonance (within 0.8 nm) yielding 
a very large negative value of a„. However, the other 
AmagicS are at far distances from the photon excitation 
energies of the resonant transitions yielding small anS. 
We also identify at least five AmagicS for the 6P3/2-651/2 
transition placed systematically between various reso¬ 
nant transitions. Similarly, several AmagicS are also lo¬ 
cated for the 65-5113/2 and 65-5115/2 transitions as seen 
in Fig. [4] in the wavelength range 300-800 nm and listed 
in Table [Tv] 

Conclusion: By determining dynamic dipole polariz¬ 
abilities of the ground nS and the first excited nP^jn^'iin 
and mD^i 2 ^^i 2 states, many magic wavelengths of the 
n5—nPi/2^3/2 and n5 —toZ 13/2^5/2 transitions of the Mg+, 
Ca’*', Sr+ and Ba’*' alkaline earth metal ions are identi¬ 
fied. Occurrence of magic wavelengths are predicted be¬ 
tween the resonant wavelengths which will be very helpful 
to perform high precision measurements in the above ions 
using both the red and blue tuned trapping techniques. 
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In a sum-over-states approach the valence correlation 
contributions can be estimated using the expression 
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where p stands for the excited intermediate states, js are 
angular momentum of the correspond state, Es are their 
energies and D is the dipole operator. It is required to 
calculate a sufficient number of atomic states so that as 
many as El matrix elements (jn||El|bp) can be evaluated 
to estimate contributions to the a.n values. Contribu¬ 
tions from the higher excited states, especially from the 
continuum, are usually small and they can be estimated 
using a lower order many-body method to a reasonable 
accuracy as suggested in our previous work [ij. We re¬ 
fer the above described dominant contribution as “Main” 
and the later contribution as “Tail” to 

Similarly, the core-valence contribution can be esti¬ 
mated by expressing as 


a. 
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where c stands for core orbitals. This basically sub¬ 
tracts the extra contributions accounted by allowing the 
Pauli’s exclusion principle of the valence electron in the 
determination of the atomic states. Typically, these con¬ 
tributions are extremely small in magnitude and can be 
estimated using a lower order many-body theory as in 
Ref. Q. 

Expressions for the core contributions is given by 
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where c and p represent sums over core and virtual or¬ 
bitals, respectively. It can be noticed that there appears 
an extra phase factor in the expression for the tensor po¬ 
larizability. Since c and p sums over a complete set of 
c and p orbitals of a closed-core, therefore core contri¬ 
bution to the tensor polarizability nullifies. Again, the 
core contributions to the scalar polarizabilities for the 
states having common core are same since its expression 
is independent of j„. 

The first four low-lying allowed transitions are included 
in the evaluations of the “Main” contributions in the 
Mg+, Ca’*' and Sr+ ions, whereas only the first three 
low-lying allowed transitions are included in the Ba^ ion. 
However, th efirst six low-lying El matrix elements are 
included in the estimation of “Main” contributions in the 
3113 / 2 , 5/2 states of the Ca+ ion. These states are evalu¬ 
ated using a linearized version of the relativistic coupled- 
cluster method with the singles and doubles ^proxima- 
tion (SD method) as described in Refs. [3, y|- Atomic 
wave function of the \jn, run) state in this method is given 

by 


|4^n) = 


1 + X Pp““p““ + 9 X Ppqabalalabtta 


pa 


pqab 


4 ^ ^ ^ PpnO’^Cln 4 " ^ ^ 

p^n pqa 
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where |d)„) is considered to be the Dirac-Eock (DF) wave 
function constructed by 

l^’n) = at |d>o) (8) 


for the DF wave function |$o) of the closed-core. In the 
above expression, at and ai represent for the creation and 
annihilation operators. The indices p, q represent for the 
virtual orbitals and the indices a, b refer to the occupied 
orbitals, ps are the amplitudes of the excitations with 
the subscripts for the corresponding orbitals. The matrix 
element of D between states Erriv) and jw'mw) states is 
evaluated in this method by J. 
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Combining El matrix elements obtained using the SD 
method with the experimental energies from the national 
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TABLE I: Calculated values of the static dipole polarizabilities of Ca^, Sr"*", Ba"*" and Ra'*" alkaline earth metal ions. These 
values are compared with the other available theoretical and experimental results. References are given in the square brackets. 
The first four low-lying allowed transitions were included in the evaluations of the “Main” contributions in the Mg’*", Ca^ and 
Sr"*" ions, whereas only the first three low-lying allowed transitions were included in the Ba"*" ion. In an exception, the first six 
low-lying El matrix elements were included in the estimation of “Main” contributions in the 81 ) 3 / 2 , 5/2 states of the Ca"*" ion. 


Mg+ 

a^°H3Si/2) 

a^°H3Pi/2) 

a^°H3P3/2) 

a‘'H3P3/2) 

a^°>(3P3/2) 

( 37 ) 3 / 2 ) 

( 37 ) 5 / 2 ) 

a^'>(3P5/2) 

Main 

34.50 

30.89 

31.18 

1.32 

185.35 

-78.41 

184.71 

-111.13 

Core 

0.48 

0.48 

0.48 

- 

0.48 

- 

0.48 

- 

Core-valence 

-1.32 X 10"^ 

0.0 

0.0 

- 

-1.12 X 10"® 

- 

-1.70 X 10"® 

- 

Tail 

6.66 X 10-2 

0.23 

0.23 

-0.16 

3.46 

-0.73 

3.46 

-1.04 

Total (Present) 

35.04 

31.60 

31.88 

1.157 

189.29 

-79.137 

188.65 

-112.18 

Others 

35.01 

31.598 

31.598 

1.163 

188.6 

-112.1 

188.6 

-112.1 

Experiment 

35.02(4) 

- 

- 

- 

- 

- 

- 

- 

Ca+ 

a^°>(45i/2) 


a‘°^(4P3/2) 

a‘"n4P3/2) 

a^°\3D3/2) 

( 37 ) 3 / 2 ) 

( 37 ) 5 / 2 ) 

a^"H3P5/2) 

Main 

72.80 

-4.97 

-2.96 

10.74 

26.23 

-16.45 

25.98 

-23.11 

Core 

3.26 

3.26 

3.26 

- 

3.26 

- 

3.26 

- 

Core-valence 

-8.85 X 10”^ 

0.0 

0.0 

- 

-7.94 X 10"® 

- 

-1.20 X 10"^ 

- 

Tail 

5.87x10"^ 

2.52 

2.52 

-0.66 

3.05 

-0.62 

2.84 

-0.82 

Total (Present) 

76.03 

0.82 

2.82 

10.08 

32.31 

-17.02 

32.05 

-23.92 

Others 

75.28 

-2.774 

-0.931 

10.12 

32.99 

-17.88 

32.81 

-25.16 

Others 

76.1(5) 

-0.75(70) 

1.02(64) 

10.31(28) 

32.0(3) 

-17.43(23) 

31.8(3) 

-24.51(29) 

Experiment [lO] 

75.3(4) 

- 

- 

- 

- 

- 

- 


Sr+ 

a™(55i/2) 

a^°H5Pi/2) 

a^°^(5P3/2) 

a‘"n5P3/2) 

a'“)(4P3/2) 

( 47 ) 3 / 2 ) 

( 47 ) 5 / 2 ) 

a^"n4P5/2) 

Main 

85.73 

-39.94 

-28.44 

11.33 

53.83 

-34.25 

53.63 

-46.33 

Core 

4.98 

4.98 

4.98 

- 

4.98 

- 

4.98 

- 

Core-valence 

-0.19 

0.0 

0.0 

- 

-1.77 X 10”^ 

- 

-2.79 X 10"^ 

- 

Tail 

1.97x10"^ 

3.70 

2.67 

-0.81 

4.95 

-1.00 

3.48 

-1.01 

Total(Present) 

90.54 

-31.27 

-20.79 

10.52 

63.74 

-35.26 

62.08 

-47.35 

Other yjj 

89.88 

-23.13 

-23.13 

- 

61.77 

- 

61.77 

- 

Other [1^ 

88.29(1.0) 

- 

- 

- 

61.43(52) 

-35.42(25) 

62.87(75) 

-48.83(30) 

Experiment [1^] 

93.3(9) 

-32.6 

-32.6 

- 

57.0 

- 

57.0 

- 

Ba+ 

a^°HGSi/2) 

a™(6Pi/2) 

a™(6P3/2) 

a^"^(6P3/2) 

a^“’(5P3/2) 

( 57 ) 3 / 2 ) 

a™ ( 57 ) 5 / 2 ) 

a^"^(5P5/2) 

Main 

114.19 

6.85 

32.00 

5.87 

39.71 

-21.92 

42.42 

-30.37 

Core 

9.35 

9.35 

9.35 

- 

9.35 

- 

9.35 

- 

Core-valence 

-0.38 

0.0 

0.0 

- 

-2.35 X 10"^ 

- 

-3.87 X 10"^ 

- 

Tail 

1.66 xlO"^ 

4.26 

4.18 

-1.17 

4.76 

-1.00 

4.80 

-1.47 

Total (Present) 

123.18 

20.46 

45.53 

4.70 

53.80 

-22.92 

56.53 

-31.83 

Other \1^ 

124.26(1.0) 

- 

- 

- 

48.81(46) 

-24.62(28) 

50.67(58) 

-30.85(31) 

Experiment [1^ 

123.88(5) 

- 

- 

- 

- 

- 

- 

- 


institute of standards and technology (NIST) database 
we determine the “Main” contributions to valence 
correlations ofonSofthe nSi/ 2 , riPi/ 2 , 3/2 and rnD 3 / 2 , 5/2 
states in the considered alkaline earth metal ions. In Ta¬ 
ble [H we list the contributions from the “Main”, core, 
valence-core and “Tail” components to the static val¬ 
ues of an for all the considered alkaline ions. The fi¬ 
nal values are compared with other available theoretical 
and experimental values. Transitions up to the 3S — 6P, 
3P — GS", 3P — 6d, 3D — 6P and 3D — 7F are included 
into the ’’Main” contribution for Mg+ ion. As seen from 


Table U the ground state dipole polarizability values ob¬ 
tained for the Mg+ ion are very close to the values esti¬ 
mated by Mitroy and Zhang Q. They evaluate the non- 
relativistic values of polarizabilities by diagonalizing the 
semi-empirical Hamiltonian in a large dimension single 
electron basis. We notice that there is a large difference 
between our polarizability results for 3 P 3/2 and 3D^/2 
states and hence the non-relativistic values calculated by 
them are unfavorable. Snow Q deduced the static po¬ 
larizability for the ground state of the Mg+ ion from the 
fine structure of high Rydberg levels. Their value 
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of 35.02(0.04) Oq agrees well with our result of 35.04 Ug. 
For Ca+ ion transitions up to the 45” — IP, AP — 7S, 
4P — 6d, 3D — 7P and 3D — 7F are included into the 
’’Main” contribution. While examining polarizabilities 
in case of the Ca’*' ion, we notice that our calculation for 
the ground state polarizability is in the agreement with 
the result obtained from experimental spectral analysis 
by Chang [I^. Also, a comparison of our calculated po¬ 
larizabilities for this ion with the calculations done using 
combination of many-body perturbation theory and SD 
scaled method by Safronova et al. @ shows very good 
agreement. Moreover, our result for the ground state is in 
consistent agreement with the results calculated by Tang 
and co-workers Q. However, a discrepancy has been ob¬ 
served between our results and other calculations among 
the P state polarizabilities this ion. We notice that the 
polarizabilities of the P states are very small due to sub¬ 
stantial cancellations between some large contributions 
to the total polarizability. The final values are, thus, 
very sensitive to these cancellations. Next, we compare 
our polarizability results for the Sr’*' ion in Table HI Tran¬ 
sitions up to the 5S — 8P, 5P — 8S, 5P — 7d, AD — 8P 
and 3D — 7F are included into the ’’Main” contribution 
for this ion. Mitroy and coworkers have used a non- 
relativistic method using a sum-over-states approach to 
determine polarizabilities of the S, P and D states of Sr+. 


It is evident from TableUthat our ground state dipole po¬ 
larizability is in agreement with their result. However, it 
may not be proper to compare their non-relativistic val¬ 
ues of the polarizabilities for the AD 3/2 and AD ^/2 states 
just by including corrections from the relativistic effects 
with the non-relativistic values agnaist our fully relativis¬ 
tic calculations. However, the static dipole polarizabil¬ 
ities calculated for the ground and AD state of the Sr+ 
ion using an ab initio relativistic coupled-cluster method 
[T^ are in close agreement with our results. There are 
no direct experimental values available to compare with 
these results. The polarizability values for the S, P and 
D states of the Sr"*" ion, derived by combining the exper¬ 
imental results with oscillator strength sums by Barklem 
and O’Mara , has a considerable discrepancy with our 
present results. In Ba’*' ion, transitions up to the 6S—8P, 
6P — 8S, 6P — 6d, 5D — 8P and 5D — 7F are included 
into the ” Main” contribution. The high-precision ground 
state polarizability measurement for Ba“'", achieved by us¬ 
ing a novel technique, based on the resonant excitation 
Stark ionization spectroscopy [T^ , is in very good agree¬ 
ment with our work. To summarize, the above analysis 
clearly justifies that our calculations of dipole polarizabil¬ 
ities in the considered alkaline earth metal ions are ac¬ 
curate enough to predict the magic wavelengths in these 
ions reliably. 
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